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1. INTRODUCTION 

Prior to launch, the space shuttle might be 
described as a very large thermos bottle containing 
substantial quantities of cryogenic fuels. Because 
thermal insulation is a critical design requirement, the 
external wall of the launch vehicle fuel tank is covered 
with an insulating foam layer. This foam is fragile and 
can be damaged by very minor impacts, such as that 
from small- to medium-size hail, which may go 
unnoticed. In May 1999, hail damage to the top of the 
External Tank (ET) of STS-96 required a rollback from 
the launch pad to the Vehicle Assembly Building 
(VAB) for repair of the insulating foam. Because of the 
potential for hail damage to the ET while exposed to 
the weather, a vigilant hail sentry system using impact 
transducers was developed as a hail damage warning 
system and to record and quantify hail events. 

The Kennedy Space Center (KSC) Hail Monitor 
System, a joint effort of the NASA and University 
Affiliated Spaceport Technology Development 
Contract (USTDC) Physics Labs, was first deployed 
for operational testing in the fall of 2006. Volunteers 
from the Community Collaborative Rain, Hail, and 



Figure 1. Hail monitor stations at Pad 39 A. 
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Snow Network (CoCoRaHS) in conjunction with 
Colorado State University were and continue to be 
active in testing duplicate hail monitor systems at 
sites in the hail prone high plains of Colorado. 

The KSC Hail Monitor System (HMS), consisting 
of three stations positioned approximately 500 ft from 
the launch pad and forming an approximate 
equilateral triangle (see Figure 1), was deployed to 
Pad 39B for support of STS-115. Two months later, 
the HMS was deployed to Pad 39A for support of 
STS-1 16. During support of STS-1 17 in late February 
2007, an unusual hail event occurred in the 
immediate vicinity of the exposed space shuttle and 
launch pad. Hail data of this event was collected by 
the HMS and analyzed. 

Support of STS-1 18 revealed another important 
application of the hail monitor system. Ground 
Instrumentation personnel check the hail monitors 
daily when a vehicle is on the launch pad, with special 
attention after any storm suspected of containing hail. 
If no hail is recorded by the HMS, the vehicle and pad 
inspection team has no need to conduct a thorough 
inspection of the vehicle immediately following a 
storm. On the afternoon of July 13, 2007, hail on the 
ground was reported by observers at the VAB, about 
three miles west of Pad 39A, as well as at several 
other locations around Kennedy Space Center. The 
HMS showed no impact detections, indicating that the 
shuttle had not been damaged by any of the 
numerous hail events which occurred that day. 

2. HAIL DISDROMETER 

Impact disdrometers (Joss, 1967) have long been 
a useful meteorological tool for measuring rainfall 
drop size distributions. The electrical impulse 
generated from a raindrop impact with the sensor 
head is converted to an estimate of drop diameter by 
means of an empirical calibration formula. The 
calibration may be a one-time procedure involving 
dropping numerous known-size-calibration drops from 
a tower of sufficient height to achieve terminal 
velocity. Alternatively, or as a supplement to the 
single drop calibration, an in situ method of comparing 
the total integrated volume of disdrometer measured 
drops per synchronized collocated tipping bucket rain 
gauge tip, as described by Lane (1997), or in the case 
of a hail disdrometer, the time integrated histograms 



3. HYDROMETEOR SIZE DISTRIBUTIONS 



Figure 2. CAD drawing of hail disdrometer 
impact sensor, showing placement of 
piezoelectric transducer. 

of hail disdrometer spectra to collocated hail pads 
(Lane, 2006) can be utilized. 

The hail disdrometer follows the design principles 
of the rainfall disdrometer. The hail sensor is 
composed of four aluminum sheet metal plates 
welded to the shape of a shallow pyramid. The 
sloped surface deflects hail, preventing multiple 
bounces from the same hail stone. A piezoelectric 
microphone element is mounted beneath and in the 
center of the pyramid (see Figure 2). The output of 
this microphone is connected to an electronic circuit 
which digitizes and processes the electrical signal, 
and then transmits a trigger pulse to one of six output 
channels. Each output channel represents a signal 
that is twice as large as the previous channel, thereby 
categorizing the hail stone into one of six sizes, from 
diameters of about 10 mm to 20 mm, in 2 mm steps. 
The six output channels are connected to six LCD 
counters, which create a permanent record of all hail 
hitting the sensor. The counters are manually reset 
after a hail storm. 



The physical connection between rainfall rate 
measured by ground instruments and rainfall 
estimated by weather radar is the drop size 
distribution (DSD). If the DSD is somehow previously 
known as a continuous function of time and space, 
radar reflectivity Z can be directly calculated from the 
spatial average of the sixth moment of the DSD. 
Using hydrometeor terminal velocity values, rainfall 
rate R can be also be calculated at any point on the 
ground if the DSD is known at ground level and air 
motion is assumed to be zero near the ground. 

Direct DSD measurements can be made at a 
point on the ground using a disdrometer. By applying 
the equations of motion for non-interacting 
hydrometeors, a spatial and temporal estimate of the 
DSD in a volume surrounding ground measuring 
stations can be obtained from a microscale array of 
disdrometers by employing a horizontal interpolation 
and vertical extrapolation algorithm such as the 3D- 
DSD (Lane, 2002a). The sampling volume is 
characterized by wind velocities, which are input 
parameters to the 3D-DSD model, composed of 
vertical (updraft/downdraft) and horizontal (advection) 
components. The 3D-DSD algorithm assumes that 
only gravity and terminal velocity due to atmospheric 
drag within the sampling volume influence 
hydrometeor dynamics, i.e., hydrometeor interactions 
are ignored. 

3.1 Gravitational Sorting of Hydrometeors 

Gravitational sorting of hydrometeors is a well 
known meteorological phenomenon. Disdrometers 
are well suited to measuring this characteristic of 
convective storms. However, the standard binning of 
disdrometer size data can mask the details of 
gravitational sorting. In order to study gravitational 
sorting it is necessary to operate a disdrometer in a 
diameter versus time (D-/) mode, e.g., McFarquhar 
(1996). 

Useful information can be extracted from D-t 
plots using a crude model of drop dynamics. For 
simplicity, the hydrometeor terminal velocity can be 
approximated by a power-law: 

v t (d)» KD p ( 1 ) 

Even though better approximations have been 
proposed, Equation (1) with p = 1/2 and K = 4.4 [m s' 1 
mm 1/2 ] is a reasonable approximation for both rain 
and hail. 

The fall time of a hydrometeor from some 
effective starting height h at t=t 0 can be expressed 
by: 

- ( 2 ) 

v t (D)-w 

where w is the vertical component of ambient air 
motion (positive direction is up). Substituting Equation 



Figure 3. Site 3 hail monitor unit at Pad 39 A. 


(1) into (2) and solving for D results in a predicted 
gravitational sorting of size diameters at the ground 
as a function of time: 


D(t) = 


/ \Up 

K(t-t 0 ) , 


( 3 ) 


where r = t-t 0 . 

Figure 4 shows a D-t scatter plot (Lane, 2002b) 
from one disdrometer of a three disdrometer array, 
located approximately 16 km north of the Melbourne 
NWS radar. This plot of D-t corresponds to a 
convective storm that occurred on the afternoon of 
July 19, 2001, beginning at 22:40 UTC. Equation (3) 
is superimposed with various values of h and w. As 
can be seen, h and w can not easily be uncoupled, 
however, the general shape of the curve is significant. 
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Figure 4. Disdrometer D-t data with calculated 
D(t) from Equation (3). 


The main point of the above discussion is to 
demonstrate that hydrometeor trajectory analysis can 
be performed, yielding useable results, based on a 
crude model that does not consider hydrometeor 
interactions. This point leads to the 3D-DSD model. 


3.2 3D-DSD Algorithm 


For the yth diameter size bin of width A D (valid for 
either rain or hail), corresponding to a hydrometeor 
diameter of Dj = jAD, the fall time A tj of the 
hydrometeor, based on its terminal velocity vj(Dj) , 
is: 

Atj = Az / v t (Dj) (4) 

where A z is the vertical distance of fall during time 
interval A tj . The corresponding horizontal distances 
A Xj and A yj , corresponding to the time interval A tj , 
are: 

AXj =u x Az!v t (Dj) (5a) 

A -i i y A z/v r (D y ) (5b) 

Together, Equations (4) and (5) describe the linear 
trajectory of a raindrop packet rjj of size Dj , 

traveling at a constant velocity , formed by the 


vector sum of the cloud advection and drop terminal 
velocities: 

Vj=u x e x +u y e y -v T (Dj)e : (6) 

where e x , e and e z are unit vectors in the x , y and 

z directions, while u x and u y are the components of 
the advection velocity. The number of drops per unit 
volume, described by the yth drop size packet 
q is related to the continuous drop size 

distribution function N(x,y,z,t 9 D ) by: 

M * 

t) * J N(x,y,z,t;D)dD ( 7 ) 

y=i o 

Using a microscale disdrometer array, 
rjj(x,y,z 0 ,t) can be estimated by interpolating the 

drop count histograms between each disdrometer of 
the array, where the interpolation is performed over x 
and y independently (z 0 is the ground level), for each 
of the yth drop size bins. The disdrometer histogram 
bins are first divided by the product of the histogram 
measurement time interval, disdrometer sensing area, 
and drop terminal velocity, in order to achieve units of 
m' 3 (i.e., drops per unit volume). Once rj.(x,y,z 0i t) 

has been estimated from the disdrometer array data, 
a interpolation to z yields the 3D-DSD model: 

rjj(x,y,z,t) * rjj(x',y',z 0 ,t') (8) 

Using Equations (4) and (5), the primed and unprimed 
coordinate systems in Equation (8) are related by: 

x'= x + u x Atj = x + u x z / v t (D) 

y'=y + u y Atj=y+u y z/v T (D) 

t'=t + Atj =/ + z/v r (D) 

According to the description established by Equations 
(8) and (9), a 3D-DSD can be estimated at any jc, y , z 
and t by transforming to a primed coordinate system 
at z 0 (ground level) where the DSD has been 
determined by measurement and interpolation. 

Reflectivity is calculated as the 6th moment of the 
3D-DSD: 

M 

Z{x,y,z,t) = Yj D< i nM,y,z,t) 

T/ (10) 

* Yj d j nj(x',y',z 0 ,n 

j - 1 

Equation (10) has strong similarities to a convolution 
sum due to the time delay (and spatial offset) caused 
by the dependence of the primed coordinates on the 
drop diameter D. The Ath disdrometer is located at 
{x k ,yk, Zo}. SO that Z(x k ,y k ,z 0 ,t) from Equation (10) 
is equivalent to the traditional reflectivity calculated 
from a point disdrometer. In order to compute a 3D- 
DSD estimate of disdrometer derived radar reflectivity 
for comparison to WSR-88D radar reflectivity, it is 
necessary to average Equation (10) over the 
appropriate extent of jc, y and z . 



Figure 5a. Melbourne WSR-88D, lowest 
three reflectivity scan elevations (top to 
bottom) at 22:09 UTC, Feb 26, 2007. 



Figure 5b. Melbourne WSR-88D lowest three 
storm relative radial velocity scans (top to 
bottom) at 22:09 UTC, Feb 26, 2007. 





Figure 6a. Melbourne WSR-88D , lowest 
three reflectivity scan elevations (top to 
bottom) at 22:14 UTC, Feb 26, 2007. 



Figure 6b. Melbourne WSR-88D lowest three 
storm relative radial velocity scans (top to 
bottom) at 22:14 UTC, Feb 26, 2007. 






4. SEVERE HAIL EVENT AT PAD 39A 

The prototype hail monitor system (see Figure 3) 
was in operation at Pad 39A at the time of the 
February 26, 2007 severe hail event. On this day the 
atmosphere over the Florida Peninsula was 
conditionally poised for scattered deep convection. 
The 1500 UTC Cape Canaveral sounding analysis 
(XMR) revealed 500 mb temperatures around -14 C 
with a mean-layer convective available potential 
energy calculated near 1500 J/kg. Despite area 
cloudiness, afternoon surface temperatures were 
anticipated to near 27 C. A southward sagging cold 
front was located over northern Florida leaving the 
KSC situated in the synoptic warm sector with vicinity 
dew point temperatures of 17 to 19 C. Opportunity for 
favorable jet dynamics aloft was evident with the right 
rear quadrant of an upper-level jet forecast to move 
across the peninsula late in the day along with an 
increasing low-level jet of 15 m/s (—30 kt) out of the 
west. Surface winds were southwest at 5 to 8 m/s (~ 
10 to 15 kt) with nearly unidirectional flow above, and 
with increasing speed shear with height. The 0-6 km 
shear was 25 to 30 m/s (~ 50 to 60 kt). The assessed 
threat was for isolated severe storms to develop late 
in the afternoon and move east at 12 to 15 m/s (~25 
to 30 kt). The associated main hazards were large 
hail, straight line winds, and excessive cloud-to- 
ground lightning. By 21:00 UTC, strong convection 
had broken out over eastern Orange County and at 
21:33 UTC a warning was issued for eastern Orange 
County and northern Brevard County for an isolated 
severe thunderstorm moving toward KSC. 
Transitioning from the rapid growth phase into the 
mature phase, an apparent "donut" signature was 
noted on the Lightning Detection and Ranging (LDAR) 
system. The LDAR system provides 3-D and 4-D 
total lightning information in proximity of the sensor 
network located on KSC and is especially useful for 
storm-cell interrogation when used in conjunction with 
weather radar. This apparent donut (e.g., lightning- 
free hole), along with an increased trend in the total 
signal count, helped discern the extreme intensity of 
the rotating updraft and bolstered forecaster 
confidence that the worst was yet to occur. 
Unfortunately, severe weather manifestation (with 
respect to large hail, stones > 19 mm at the surface) 
occurred before the storm could move offshore and 
away from the shuttle launch pad and KSC 
compound. 

4.1 Melbourne WSR-88D 

Volumetrically, a once elevated core of very high 
reflectivity (~70 dBZ) rapidly descended to the surface 
as shown in the lowest elevation slice at 22:09 UTC of 
Figure 5a. Large hail aloft was no longer being 
suspended by the intense updraft, signaling a shift to 
the downdraft dominated phase of this severe local 
storm. This phase represents the surface impact 
phase; the dangerous phase in this scenario. The 
storm-relative velocity (Figure 5b) shows broad 


cyclonic vorticity surrounding the main core and 
denoting some descending angular momentum 
influenced from organized mid-level rotation noted 
several volume scans earlier (the organized mid-level 
rotation and storm tilt with height allowed hail 
embryos to be suspended for longer periods with 
preferred trajectories for greater growth opportunity). 
At 22:09 UTC, although radial magnitudes were 
notably high, the area was broad (not an intense gate- 
to-gate couplet; not a significant tornado threat). A 
flanking gust front moving east southeast across the 
Pad 39A was situated near and to the immediate 
southwest, with subsequent low-level convergence 
just ahead of it (promoting additional cellular growth). 
Relative to Pad 39A, the strongest winds and largest 
hail appeared to hit between 22:09 UTC and 22:14 
UTC (Figure 6). Accelerated winds of 25 to 30 m/s 
(—49 to 58 kt) were noted south of the pad, as well as 
a turbulent eddy (or possible gustnado/funnel). 

4.2 Hail Pad Data 

Data from the Pad 39A hail monitor system 
recorded hail diameters up to 20 [mm]. Collocated hail 
pads (Styrofoam covered with aluminum foil, based 
on the CoCoRaHS design and fabrication methods) 
suggest hail sizes up to 35 [mm] may have occurred 
at Sites-2 and -3 (see Figures 8 and 9). A hail pad 
calibration method developed during the 1978 Alberta 


Figure 7. Site-1 hail pad after Feb 26 hail event ; 
(bottom) foil removed, with lxl [in] grid. 





Figure 8. Site-2 hail pad after Feb 26 hail event; 
(bottom) foil removed, with lxl [in] grid. 


Hail Project relates the dent diameter D to the hail 
diameter D'as an empirical second-order polynomial 
(Strong, 1977): 

D' = a 0 +a t D + a 2 D 2 (11) 

where a 0 = 0.38 [cm], a x - 1.11, and a 2 = -0.04 
[cm -1 ]. Other calibration formulas have since been 
developed (Lozowski, 1978), as well as calibrations 
for hail pads made from alternate materials, such as 
foam rubber covered with aluminum foil (Waldvogel, 
1978b), and Styrofoam covered with white latex paint 
(Long, 1980). 

Two of three hail pads at Sites-2 and -3, 
experienced hail driven by violent horizontal winds 
thus shredding the foil (see Figures 7 through 9). A 
small number of hail pads with shredded foil are seen 
every year by the Colorado State University managed 
CoCoRaHS network of hail observers, usually only 
occurring when numerous stones of moderate to large 
size are driven by strong winds. Out of the many 
powerful hail storms that occur in Colorado each year, 
only about one out of every 100 or so hail pads are 
retrieved with significantly torn foil. 

The large impact angle deviation from normal 
incidence makes the 30-35 [mm] hail pad deduced 
maximum diameters at Sites-2 and -3, most likely an 
overestimate due to the gouging effect as opposed to 
the denting mechanism of a normal incidence impact 



Figure 9. Site-3 hail pad after Feb 26 hail event; 
(bottom) foil removed, with lxl [in] grid. 


at terminal velocity. However, the fact that gouging 
occurred, is relevant to the total kinetic energy and 



[kt] peak from 270°, with a duration of about one 
minute. 

hail damage potential on the shuttle ET foam. 

4.3 Hailstone Spectra 

Rainfall disdrometer data leads directly to an 
estimation of the corresponding rainfall DSD. Rainfall 


DSD estimates can then be used to compute an 
equivalent radar reflectivity via a discrete evaluation 
of the 6 th moment of the DSD. An estimation of 
rainfall rate can also be performed via discrete 
integration (summation of histogram bin products) 
using raindrop terminal velocity estimates. 


Table 1. Hail disdrometer size channel counts. 



HM Site-1 

HM Site-2 

HM Site-3 

Chan, 

A* 

Cnt, 

D k 

Cnt, 

A* 

Cnt, 

k 

[mm] 

h k 

[mm] 

h k 

[mm] 

h k 

1 

10.7 

28 

11.4 

124 

9.6 

103 

2 

12.8 

45 

13.9 

39 

12.0 

89 

3 

14.5 

11 

15.9 

8 

13.9 

24 

4 

16.2 

0 

17.9 

2 

15.8 

5 

5 

17.9 

0 

19.9 

0 

17.7 

1 

6 

>18.7 

0 

>20.9 

0 

>18.6 

0 


Table 2. Hail disdrometer counts per 60 [s] interval. 



HM Site-1 

HM Site-2 

HM Site-3 

Logger 
clock offset 
correction 
= 0[s] 

Logger 
clock offset 
correction 
= -24 [s] 

Logger 
clock offset 
correction 
= +6 [s] 

Time 

Tot Cnt, 

Tot Cnt, 

Tot Cnt, 

‘j 

Wj 

Wj 

Wj 

22:09 

0 

0 

0 

22:10 

0 

0 

1 

22:11 

4 

38 

96 

22:12 

57 

112 

109 

22:13 

18 

10 

6 

22:14 

3 

10 

7 

22:15 

1 

2 

0 

22:16 

0 

1 

3 

22:17 

1 

0 

0 

22:18 

0 

2 

1 

22:19 

0 

0 

0 


In the case of hailstone spectrometers (see for 
example Waldvogel, 1978a), a hail impact 
disdrometer, such as that used in the Pad 39A hail 
monitor system is a subset of that class of 
instruments. Hail size distributions (HSD) can also be 
computed from hailstone spectra, along with an 
equivalent radar reflectivity (of hail only). Instead of a 
hail fall rate equivalent to rainfall rate, a hailstone 
kinetic energy flux is more often the quantity of 
interest, since the kinetic energy flux or total 
integrated energy is more relevant to damage caused 
by hail, e.g., Visser (2000). 

Kinetic energy flux E„ [J m' 2 s' 1 ] is defined by 
Waldvogel (1978a) as: 

E H = np . \N H {D)D\\ dD < 12 ) 
* 12xl0 6 j " T 


where p [gm cm' 3 ] is the density of hail, AWD) [m -3 
mm' 1 ] is the hail size distribution (HSD), and v T is hail 
stone terminal velocity, such as that approximated by 
Equation (1). The total kinetic energy E H [J m' 2 ] per 
unit area is just the time integral of Equation (12) over 
the duration of the hail event. 

The kinetic energy associated with the Feb 26 
hail event at Pad 39A can be computed from the hail 
disdrometer data shown in Table 1 and by using a 
discrete version of the time integral of Equation (12): 


e h = 


np 1 

12x10 6 T 


\KD\v 


(13) 


where h k is the &th channel count corresponding to D k 
from Table 1, and v Tk = v T (D k ) from Equation (1). 

A s is the disdrometer sensor area, equal to one sq ft 
or 0.093 [m 2 ]. 

The result of evaluating Equation (13) for each of 
the three hail monitor sites is summarized in Table 3. 
These values are reasonable when compared to 
published kinetic energy flux, e.g. Visser (2000), of 
E h « 1 8 [J m' 2 s' 1 ] for radar reflectivity values of 66 
dBZ. The energy flux is estimated in Table 3 by 
dividing the total energy density by A/ = 120 [s], which 
is the approximate time interval of the majority of the 



Figure 11. Hail size distribution for the Feb 26, 2007 
event calculated using Equation (13) with Pad 39 A 
hail data from Table 1. Note comparison to HSD 
models from Federer (1975) and from Auer (1972). 


hail fall. Since the hail disdrometer at Site-1 was 
located on the leeward side of the pad and was 
therefore partially sheltered from the hail event, it is 
not surprising to see an energy flux value that is 
smaller than the other two sites. 


Table 3. Total kinetic energy and estimated flux of 
Feb 26 hail event, with At = 120 [sj. 



HM Site-1 

HM Site-2 

HM Site-3 

Eh [J 

nr 2 ] 

99.6 

209.5 

194.4 

E h « 

Ejl 

At 

0.83 

1.75 

1.62 

[J m' 2 

s 1 ] 





The hail size distribution, Nh{D) [m' 3 mm' 1 ], can be 
calculated from the histogram data in Table 1 and the 
following equation: 


N H (D k ) 


K 

A s v Tk At AD 


(13) 


where h k is again the &h channel hail impact count 
corresponding to D k in Table 1; v Tk =v T (D k ) from 

Equation (1); A s = 0.093 [m 2 ] ; At = 240 [s] (a longer 
event interval than the energy flux case); and AD * 
2.0 [mm]. The result is plotted in Figure 11 and for 
comparison, Marshall-Palmer model plots of rainfall 
DSDs for various values of rainfall rate, as well as 
some published hail size distribution (HSD) models. 

The HSD plots from Pad 39A in Figure 11 agree 
reasonably well with Federer (1975). Agreement with 
Auer (1972) below hail diameters of 15 [mm] is also 
fair. The droop and cut-off of the Pad 39A HSD below 
12 [mm] is most likely due to the calibration of the hail 
disdrometers and software cut-off limit which was set 
to aggressively discriminate rain from hail, i.e., 
adjusted for a low probability of false detection. This 
threshold may be moved to a lower value with some 
increased risk of false detections from rain. 


4.4 3D-DSD Plots 

Using the hail disdrometer array data from Tables 
1 and 2, as well as Equation (10) as a starting point, a 
disdrometer derived radar reflectivity can be 
computed for comparison to the NWS WSR-88D (see 
Figure 12). Several steps need to be taken in order to 
optimize this comparison: 

• Vertical (updraft/downdraft) and horizontal 
(advection) wind velocities are used to refine the 
calculations of Equation (10), as discussed by 
Lane (2002b). 

• The time and location of the corresponding WSR- 
88D reflectivity bin center is used in Equation 
(10) for each computation of x, y, z, and t in the 
disdrometer 3D-DSD plot of equivalent 
reflectivity. 


Melbourne NEXRAD Reflectivity Hail Monitor - 3D-DSD 



Melbourne NEXRAD Reflectivity Hail Monitor - 3D-DSD 



Melbourne NEXRAD Reflectivity Hail Monitor - 3D-DSD 



Melbourne NEXRAD Reflectivity Hail Monitor - 3D-DSD 
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Melbourne NEXRAD Reflectivity Hail Monitor - 3D-DSO 



Figure 12. Comparisons of WSR-88D radar 
reflectivity and equivalent radar reflectivity from 3D- 
DSD algorithm calculated using Equation (10) and 
hail disdrometer data from Tables 1 and 2. 
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Melbourne NEXRAD Reflectivity Hail Monitor - 3D-DSD 
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Figure 13. Same comparisons as 3 rd and 4 th set of 
Figure 12, but with the addition of a rainfall 
background based on the Marshall- Palmer drop size 
distribution; (top) background R = 130 [mm K 1 ]; 
(bottom) background R = 10 [mm K 1 ]. 


• Equation (10) is computed for numerous adjacent 
values of x, y, and z, then averaged in order to 
approximately match radar reflectivity averaged 
over a radar bin. 

• A Marshall-Palmer rainfall DSD background is 
added to the argument of the summation in 
Equation (10), as shown in Figure 13. 

The plots of Figure 12 and 13 are centered over 
Pad 39A, and represent a 1 [km] x 1 [km] horizontal 
extent with a 0.5 [km] vertical extent. Figure 12 
shows comparisons of WSR-88D radar reflectivity and 
equivalent radar reflectivity from the 3D-DSD 
algorithm of Equation (10), where the plotted 
disdrometer derived radar reflectivity is due to the 
contribution from hail only. 

Figure 13 shows the same comparisons as the 
22:14 UTC and 22:19 UTC scans of Figure 12, but 
with the addition of a rainfall background based on the 
Marshall-Palmer drop size distribution. The 22:14 
UTC scan uses a background R = 130 [mm h' 1 ]. The 
22:19 UTC scan uses a background R = 10 [mm h' 1 ]. 

5. CONCLUSIONS 

The hail swath that damaged STS-117 on Pad 39A, 
just brushed the pad structures. This was very 
fortunate, since the main core of the hail was roughly 
1 [km] or so to the south. According to data analysis, 
the worst part of the hail swath which contained the 


larger hail, ended just beyond the pad structures, and 
was composed of only small hail stones at the north 
perimeter fence line, about 0.5 [km] from the pad 
center. The unfortunate part was that the hail was 
driven by a 60 [kt] horizontal wind for roughly 1 [min]. 
Those conditions undoubtedly led to the hail damage 
on the vehicle’s ET. Based on radar and hail 
disdrometer data, it is clear that the hail at pad 39A 
was embedded in intense rainfall of a high rainfall 
rate. 

The HMS consisting of three hail disdrometers, 
has proven to be a useful hail monitoring and post 
analysis tool for shuttle launch pad ground support. 
The 3D-DSD algorithm, coupled with other instrument 
data, such as Kennedy Space Center wind tower and 
Melbourne NWS radar, provide a high resolution 
picture of hail events on a spatial resolution of meters, 
and a temporal scale of seconds. 
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